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We review the status of particle interferometry in ultra-relativistic 
nucleus-nucleus collisions. The theoretical focus is on the model- 
independent space-time interpretation of HBT radius parameters and 
its extension to the geometrical and dynamical asymmetries generated 
in finite impact parameter collisions. On the experimental side, we give a 
complete account of all presently available data for beam energies above 
2 AGeV. We discuss what these data imply for the dynamics and space- 
time extension of the collision region and for the condition under which 
particles decouple from this region. 
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1. Introduction 

One aspect in which a theory of high-energy nucleus- nucleus collisions has to 
go beyond conventional nuclear and high-energy physics is in understanding 
the space-time extension of the collision region and its collective dynamical 
evolution. On the microscopic level, basic quantities governing the reaction 
dynamics such as the probability that a particle can escape the collision 
region without further interaction ("freeze-out") or that it participates in 
three-body interactions, clearly depend on the space-time extension of the 
collision region and the relative velocity of the particle with respect to the 
medium. On a macroscopic level, the study of the equation of state of hot 
and dense matter rests on assessing the energy density and pressure attained 
in nucleus- nucleus collisions. The operational definition of these quantities 
clearly involves the measurement of volumes and the determination of the 
degree of collectivity of the expansion. 

It has been long recognised in both astrophysicsEJH and particle 
physicscZl that correlations between identical bosons give access to the ge- 
ometry of particle emitting sources. While statistical requirements for par- 
ticle correlation measurements are significant, interferometric techniques 
are often the method of choice when conventional space-time measurement 
techniques (like paralaxe measurements or scattering with external probes) 
are not feasible. In the context of nucleus-nucleus collisions, the use of 
identical two-particle "Hanbury-Brown Twiss" (HBT) correlations was first 
forcefully advocated in the_L27Qs.in a series of oiDpfiering works published 
independently in the EastEI'E3'E3 and the WestEil'Ell. Those were followed 
by the first detailed discussions of how to separate quantum interference 
effects from final state interactions e3H. By now, particle interferometry 
provides the most direct and most detailed space-time picture of the late 
stage of nucleus-nucleus collisions. 

2 AGeV, the lowest beam energy delivered by the Alternating Gradi- 
ent Synchrotron (AGS) at the Brookhaven National Laboratory (BNL), is 
a natural choice for the lowest energy included in a review on relativistic 
nucleus- nucleus collisions. Non- relativistic collisions at lower beam energies 
(up to several 100 AMeV) create evaporative sources with long lifetimes 
('^ several 100 fm/c) which are dominated by proton and neutron emission. 
Strong final state effects dominate correlation measurements. A comprehen- 
sive review of the rather different physics at these lower energies existal3. 
Above 2 AGeV, pion abundances become important, source lifetimes are 



4 



B. Tomdsik and U.A. Wiedemann 



much shorter (~ 10 fm/c) and quantum interference effects dominate over 
the strong (but not necessarily over the electro-magnetic) final state interac- 
tions. Moreover, the increase of pion multiplicity per event with increasing 
beam energy helps to overcome the major obstacle of particle interferometry 
at lower energies, namely the lack of statistics. 

The increase in event multiplicity and the accumulation of increasing 
event samples lead throughout the 1990s to a rapid refinement of analysis 
tools. Early analyses often show "one-dimensional" two-particle correla- 
tors as functions of one relative momentum variable only. These were soon 
superseded by "three-dimensional" representations. Then, the increasing 
statistics allowed to resolve "three-dimensional" correlation functions with 
respect to the average transverse pair momentum and rapidity. Only re- 
cently, the next step of these refinements was taken when the first measure- 
ments of HBT with respect to the orientation of the reaction plane became 
available. Parallelling this experimental progress and sometimes anticipat- 
ing it was a theoretical effort which related systematically these more and 
more differential interferometric measurements to more and more detailed 
geometrical and dynamical properties of the particle emitting source,_&tJ 
end of the CERN SPS heavy ion program in 2000, several reviewst2j'Li 
gave detailed summaries of these developments. 

We were asked to contribute to the present review volume only two 
years later, at a time when HBT measurements from RHIC are published 
but no thorough theoretical analysis of these findings is yet available. In this 
situation, we decided to emphasise in a topical review those two aspects of 
particle interferometry which are likely to play an important role in the fur- 
ther discussion of RHIC data and which are not yet adequately represented 
in reviews. On the experimental side, our focus will be on a comprehensive 
overview of all data from AGS, SPS and RHIC. This is timely not only 
because several experiments at the AGS and SPS just finalised their HBT 
analyses. It is also of obvious use for a discussion of the energy-dependence 
of HBT correlation measurements which plays an important role in the 
HBT studies at RHIC. On the theoretical side, our presentation reflects 
the fact that the first measurements of azimuthally dependent HBT ra- 
dius parameters in finite impact parameter collisions at AGS and RHIC 
became available after 1999. That these novel observables may provide a 
complementary space-time picture to the large elliptic flow measured at 
RHIC gives an additional motivation for discussing their theoretical basis 
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in detail. 

This review is organised as follows: In section |[ we introduce the formal- 
ism for calculating identical two-particle correlations from particle emission 
functions. In particular, we discuss in detail the case of non-central colli- 
sions. In section we briefly summarise models for the particle emitting 
source which have been used in calculating HBT correlations from ultra- 
relativistic heavy- ion collisions. Finally, in section ^, we give a complete 
account of the experimental situation and we discuss shortly the most im- 
portant features of the presently available data. Since most data are for 
central collisions, readers primarily interested in the data summary may 
jump from the introductory sections 2.1 and 2.2 directly to section^. The 
extension to non-central collisions in section |2.3| and the corresponding ex- 



perimental overview in section 4.5.5 can be read independently. Also, the 
short overview of model calculations in section |^ allows for an independent 
reading. 



2. The Gaussian source formalism 

The two-particle momentum correlator between identical bosons with mo- 
menta pi and p2 is defined as the quotient of two-particle and one-particle 
spectra^. 



Cil,K)^^S^. (1) 

Experimentally, this ratio is constructed for event samples, normalising ac- 
tual pairs from the same event by mixed pairs from different events. Usually, 
the two-particle correlator is written in terms of the corresponding relative 
and average pair momenta 

q=pi-p2, (2) 
K^l{pi+P2). (3) 

If the particle spectra entering (|l|) are corrected for final state interactions, 
then the twO|4iai.ti£le.^CQp£lator C{q, K) is related to the emission func- 
tion S{x, K) tMmmml which is the Wi gner phase-space density of the 



Throughout this review, we use itaUcs for four-vectors and bold for three-vectors. 
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particle emitting system 

\JSxS{x,K)\^ 

For a derivation of in different formalisms (coherent source formalism 
and Gaussian wave-packet formalism), we refer to the reviewt22l. Equa- 
tion (^) is approximate since it is based on the smoothness approxima- 
tion S{x,K — ^q)S{y,K + ^q) « S{x,K)S{y,K) valid for small rela- 
tive momenta. Moreover, the following discussion of will regularly em- 
ploy the on-shell approximation Kq « y/K"^ + m? which neglects the off- 
shell components of K. Corrections to these approximations were stud- 
ied systematicallyllj and were found to be small except for especially con- 
structed toy models whose phase-space volumes were very small (i.e. com- 
parable to the lower bound coming from the Heisenberg uncertainty rela- 
tion). For particle correlations in nucleus- nucleus collisions, equation (^ is 
a well-defined starting point. The emission function S{x, K) can be viewed 
as the probability that a particle with momentum K is emitted from the 
space-time point x in the collision region. 

The aim of HBT two-particle interferometry is to extract from the exper- 
imentally measured correlator C(<7, K) as much information about S{x, K) 
as possible. This includes i) information about the geometrical extent of 
the collision region at the time of last hadronic scattering (freeze-out), ii) 
information about the freeze-out time and emission duration of the source, 
iii) information about the collective velocity of the expanding collision re- 
gion in the directions parallel and orthogonal to the beam, iv) information 
about the azimuthal dependence of the geometrical and dynamical prop- 
erties of the collision region in finite impact parameter collisions. From 
these primary measurements, various derived quantities can be extracted, 
in particular the particle phase-space density at freeze-out. 

2.1. Gaussian correlators in terms of space-time variances 

Experimental data of two-particle correlations are usually parametrised by 
a Gaussian ansatz in terms of HBT radius parameters Rij{K) and the 
A-intercept parameter 



C{q,K) = 1 + \{K) exp 



(5) 
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Here, the indices i, j run over three of the four components of q. The 
fourth component of q in (|^) is fixed by the requirement that the final state 
particles are on-shell 

q-K^l{pl-pl)^O^qo = ^. (6) 

Different choices of the three independent components of (7 = [qg , q) corre- 
spond to different Gaussian parametrisations. 

The connection between the HBT radius parameters and the space- 
time structure of the source is based on a Gaussian approximation to the 
emission function 

'2 



S{x,K) « S{x{K),K)exp 



-i^^{K)B^,{K)5:-{K) 



(7) 



The space-time coordinates x^i in are defined relative taJhe "effective 
source centre" x{K) for bosons emitted with momentum /dlaa 

i^'{K)^x^'~ x" (K) , x" [K) = (a;^) [K] , (8) 

where (...) denotes an average with the emission function S{x, K): 

, , fd^xf(x)S(x,K) 

The choice 

{B-')^,{K) = {i^S:,){K) (10) 

ensures that the Gaussian ansatz (|^) has the same rms widths in space-time 
as the full emission function. Inserting into the basic relation (^), the 
correlator takes a Gaussian form, 

C{q,K) = l + cxp[-q^q,{i^i-'){K)] . (11) 

Since the correlator depends only on the relative distances i'^ with respect 
to the source centre, no information can be obtained about the position 
x{K) of the centre of the emission. 

Finally, we comment on X{K) in eq. which parametrises the inter- 
cept of the correlation function for g = 0. This A-parameter is unity for a 
chaotic source and it is smaller than 1 for a source with partially coherent 
particle emission. In practice, however, there are many other reasons why 
deviations from C(a~ 0) = 2 may be observed. For example, contributions 
to pion productiorO from long-lived resonances lead to an extended source 
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Fig. 1. The out-side-longitudinal osl coordinate system. 



component which is reflected in a narrow peak of the correlation function. If 
this peak is narrower than the experimental resolution in relative momen- 
tum, a reduced intercept is observed. Also, misidentified particles in the 
sample of identical particle pairs reduce the correlation strength A(K). In 
addition to this, many technical details of the analysis - such as Coulomb 
correction or finite momentum resolution - can influence the value of A. 

2.2. Central (azimuthally symmetric) collisions 
2.2.1. The Cartesian Bertsch-Pratt parametrisation 

In the out- side-longitudinal (osl) coordinate system, defined in Fig. |l|, the 
longitudinal {long) direction points along the beam axis. In the transverse 
plane, the oui direction is chosen parallel to the transverse component of the 
pair momentum the remaining Cartesian component denotes the side 
direction. In this way, the out-axis is rotated independently for each particle 
pair. The azimuthal symmetry of the particle source in central collisions 
then translates in the osl coordinate system into a reflection symmetry 
with respect to the side-direction: 



S{x,y,z,t;K±,Ki) = S{x,-y, z,t; K±, Ki) 



(12) 
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This y — !■ —y symmetry translates to a — symmetry of 



the two-particle_jaj)xijelation function (11). The Cartesian Bertsch-Pratt 
parametrisationES'EHla exploits this reflection symmetry to write the two- 
particle correlator as 

C{q, K) = cM-Rl{K) ql - Rl{K) ql - Rf{K) qf - 2 RUK) q^qi] . (13) 

This parametrisation is based on the three Cartesian spatial components 
qo (out), qs (side), qi {long) of the relative momentum q. The temporal 
component is eliminated via the mass-shell constraint (||) 

q"^l3-q, /3-(/3i,O,/30 , (14) 

where (3 = K/Kq. To express the HBT radius parameters Rlf{K) in 
terms of space-time variances {x^x^){K)^ we insert the on-shell constraint 
= (3 q into the correlator (^l]) , and compare the result with the Gaussian 
parametrisation ^ of the correlator. One then finds that the HBT radius 
parameters measure different combinations of the spatial and temporal ex- 
tent of the collision systemEScj: 

Rl{K) = {f){K), (15) 

Rl{K) = {{i~(3jf){K), (16) 

RKK)^{{i-Pm{K). (17) 

Rli{K) = {{i - (3j){S - m){K) , (18) 

RlAK) = 0, (19) 

RUK) = Q. (20) 

The y ^ —y reflection symmetry of the emission function implies that the 
three space-time variances {x^Xi,){K) linear in y vanish. For the azimuthally 
symmetric case, there are only seven non-vanishing space-time variances 
which combine to four non- vanishing HBT-radius parameters R^^ (K) . This 
information is tabulated in the first column of Table ^. 

In general, symmetries of the emission function translate into constraints 
for the space-time variances and thus lead to simpler expressions for the 
HBT radius parameters. Of particular relevance is the limit K± 0. In 
this limit, no transverse vector allows to distinguish between the out- and 
side-components, and space-time variances thus become invariant under the 
exchange x ^ y. This implies that {x^){K± = 0) = {y^){K± = 0). Also, 
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the off-diagonal terms {zx){K± ~ 0) and {tx){K± = 0) vanish. As a conse- 
quenee, the HBT radius parameters show the following limiting behaviour 

lim Rl(K)= lim RliK) , (21) 

lim RlAK) = 0. (22) 

Further symmetries are sometimes helpful to gain physical intuition. For ex- 
ample, longitudinal boost-invariance implies a z — > — z reflection symmetry 
of the emission function which is approximately satisfied near mid-rapidity. 
This implies that the space-time variances linear in z vanish. In the longi- 
tudinally comoving system (LCMS) 

LCMS: l3i = Q, (23) 

this leads to the further simphfications R^{K) = {P){K) and R^i{K) = 0. 



2.2.2. Interpretation of space-time variances 

The space-time variances {Xfj_x^){K) depend on the pair momentum K. 
To understand their physical meaning, consider an observer who views a 
strongly expanding collision region. Some parts of the collision region move 
towards the observer and the particle spectrum emitted from those parts 
will appear blue-shifted. Other parts move away from the observer and ap- 
pear red-shifted. Thus, if the observer looks at the collision system with a 
wavelength filter of some frequency, he sees only part of the collision re- 
gion. Adopting a notion coined by Sinyukov, the observer sees a "region of 
homogeneity" . In HBT interferometry, the role of the wavelength filter is 
played by the pair momentum K. The direction of K corresponds to the 
direction from which the collision region is viewed. The modulus of K char- 
acterises the central velocity (3 = K / Kq of that part of the collision region 
which is seen through the wavelength filter. This situation is illustrated in 
Fig. 1^. Depending on the direction and modulus of different parts of the 
collision region are measured. The centre of a region of homogeneity (x^J 
depends on K and lies typically between the center of the collision region 
and the observer. In the Gaussian approximation, the region of homogene- 
ity is described by a four-dimensional space-time ellipsoid centered around 
{Xfj,){K) and characterised by S{x, K) in eq. (0). The widths of this region 
of homogeneity correspond to the space-time variances {xi^iX^){K), see eq. 
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flow determines 
and size of 




Fig. 2. Sketch of the particle emission region for a colhsion with collective dynamical 
expansion. The size and orientation of the average pair momentum K serves as wave- 
length filter: depending on K (indicated by arrows), HBT measurements are sensitive 
to different parts of the collision region. 



1C|). Thus, HBT radius parameters give access to the space-time variances 
{x^Xu){K) but they do not depend on the "effective source center" {x^). 



2.2.3. The Yano-Koonin-Podgoretskii parametrisation 

The Yano-Koonin-Podgoretsfcii parametrisation ehminates the mass-sheii 
constraint K -q — Q in (^) by choosing the reiative momentum components 
q± = q"^ + q ^, rp ax \c\ qi as independent. The corresponding Gaussian 



ansatz reads" 



Ciq, K) = l + X{K) exp -Rl{K) ql - Rj{K){qf - {q"^ 



(24) 



Here, the 4-veiocity U{K) has oniy one fC-dependent iongitudinal spatial 
component, 

1 



U{K) = 



(l,0,0,i;(K)) 



(25) 



The combinations of relative momenta {qf — (q°)^), {q-U{K))'^ and q^ ap- 
pearing in (^) are scalars under longitudinal boosts. Therefore, the three 
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YKP fit parameters R\{K), Rq{K), and are longitudinally boost- 

invariant and do not depend on the longitudinal velocity of the measure- 
ment frame. The fourth YKP parameter is the Yano-Koonin (YK) velocity 
v{K). The corresponding rapidity 

v{Ky 



1 



■In 



1 



(26) 



2'" \l-v{K) 

transforms additively under longitudinal boosts. 

Identifying the combinations of relative momenta (qf — (g") ), {q-U{K))'^ 
and in the Gaussian correlator (^l|), one can relate the YKP parame- 
ters to space-time variances in complete analogy to th e Ca rtesian case. Ex- 
plicit expressions can be found in the original papersE3'E£il. Since the ansatz 
( |2^ ) uses four Gaussian parameters, it provides a complete parametrisation 
for azimuthally symmetric collisions which is equivalent to the Cartesian 
Bcrtsch-Pratt parametrisation. The Gaussian HBT radius parameters of 
the Caitesian parametrisation can be related to the YKP fit parameters 



d2 d2 



r)2 d2 r?2 



1 



i?f = (1 - f3f) + 



\2 



v^Rl 



(27) 
(28) 



(A 



Rn + R 



(29) 
(30) 



l-v'^{K) 

Rl, = /3x {-mi + /f'J(^) [rI + Rl 

The s pace- t ime interpretation of YKP parameters is studied in 
detail! 



1 The consistency check (|27])-(30) of measured YKP and Carte- 
sian parameters is mandatory for a space-time interpretation of YKP pa- 
rameters since known technical prohlsms with the YKP parametrization 
can affect fit parameters significantlyE^. 



2.3. Collisions at finite impact parameter 
2.3.1. Choice of coordinate system 

For central collisions, azimuthal symmetry of the collision region allows to 
use the osZ-coordinate system which is oriented differently for each particle 
pair, such that its out-axis is parallel to the average pair momentum. 



HBT from AGS to RHIC 



13 



reaction 



plane 




out 



X 



► b 



Fig. 3. The impact parameter fixed coordinate system is obtained by rotating the out- 
side-longitudinal coordinate system by the azimuthal angle $ = A{Kj_,b). 

At finite impact parameter, HBT radius parameters Rfj{K) depend 
not only on the modulus of the transverse pair momentum \K\, but also 
on its orientation in the transverse plane which is defined with respect to 
some pair-independent direction in the laboratory systemlloS. In prac- 
tice, one introduces the impact parameter fixed cooTdinede system, for which 
the transverse coordinate x lies within the reaction plane; the transverse 
coordinate y lies orthogonal to it, see Fig. ^ The out-side-longitudinal and 
the impact parameter fixed system are related by a rotation with respect 
to the angle $ = Z{Kj_,b) 



For non-central collisions, the emission function is written in the impact 
parameter fixed coordinate system. 

2.3.2. HBT radius parameters from space-time variances 

In eq. ^ , the emission function entering the two-particle correlator is writ- 
ten in the osZ-system. For an emission function in the impact parameter 




such that after rotating to the osZ-system we have 




(32) 
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fixed system, the correlator reads: 

C{K, g) = 1 + I (e"J-(^-»'-^-/3*)) I . (33) 
From this, the HBT radius parameters can be calculatedEl: 



K{K) = - 



dq^ dqj 



q=0 



= {[{V^A)^ - iV^f3)^i][{V^x), - {V^!3)jt]) . (34) 

As for the azimuthally symmetric caae, the HBT radius parameters show 
imphcit and explicit X-dependencea^S: 

Rl{K±, $, Y) = (i^) sin^ $ + {f) cos^ $ - {iy) sin2$ , (35) 
Rl[K^,<^,Y) ^ (£2)cos2$+ (y2)sin2$ + /32^P) 

-2(3^{tx) cos$ - 2l3±{iy) sin$ + (xy) sin 2$ , (36) 
Rl,{K^,<P,Y) = (xy) cos2$ + isin2$((y2) _ (£2^) 

+/3± (iS;) sin $ - /3± (iy) cos $ , (37) 

Rf{K^,^,Y) = {Cz-f3iif), (38) 

i?^i(ifi,$,y) = ((z-Ai)(Scos$ + ysin$-/3it)), (39) 

i?2;(if^,$,y) = ((f-Af)(ycos$-isin$)). (40) 

On the r.h.s. of these equations we only displayed the explicit ^-dependences 
which are a geometrical consequence of rotating the owt-axis from the di- 
rection of the impact parameter b to the direction of the average pair mo- 
mentum K±. The implicit ^-dependence of the space-time variances, 

(i^i^) = {Xf,x,,){K±, $, Y) , (41) 

characterises the dynamical correlations between the size of the effective 
emission region and the azimuthal direction in which particles are emitted. 
The role of explicit and implicit contributions is illustrated in Fig. ^. 

2.3.3. Symmetries satisfied by space-time variances 

The implicit ^-dependence of space-time variances can be characterised in 
terms of the Fourier coefhcients {■■■)%, and {■ ■ ■)n 

oo 

(5^5^)(4>) = ^(i^i^)^, cos(n$) -f- (i^ii.)^ sin(n$) . (42) 

n=0 
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low K, 




(]) = - 71/2 













high K, 















Fig. 4. Schematic picture of the role of exphcit and imphcit t/i-dependence in the trans- 
verse plane. At small transverse pair momentum, the orientation of the region of ho- 
mogeneity follows the global geometry, and the explicit ^-dependence is expected to 
dominate. At large transverse pair momentum, the orientation and shape of the region 
of homogeneity can deviate significantly from the global geometry. In this case, the im- 
plicit ^-dependence of HBT radius parameters becomes significant, see text. 



Several symmetries allow to constrain this most general implicit <i>- 
dependenceS'0ii: 



(1) Mirror symmetry with respect to the reaction plane: 

This symmetry of the emission function holds at both mid and forward 
rapidity. It implies that space-time variances linear in y change sign 
under ^ — while space-time variances which are not linear in y 
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are invariant under the mirror symmetry, 

{Xfj_x^){^) ^ ~{x^^Xy){~<^) , for either = y or 5^ = 2/ , (43) 
(i^i^) ($) = {x^,x^) (-$) , for 5^ , ^ y , (44) 

mm^mi-^). (45) 

(2) Point symmetry at mid rapidity: 

For coUisions between equal mass nuclei, the total emission region is 
point symmetric under (i, y, z) —f (— i, —y, —z) at mid-rapidity. This 
is a symmetry under $ — ^ $ + tt. At mid-rapidity, we thus have the 
additional requirement 

{x^J,Xy){^') = —{x^J,Xy){^ + vr) , for either x^^ — t or Xy —i, (46) 
>(*) = ( )($ + 7r), for i^.iy^t, (47) 
(?)($) = (P)($ + 7r). (48) 

At forward rapidity, this point symmetry does not constrain the az- 
imuthal dependence since the same symmetry maps Ki — > —Ki on 
which the space-time variances depend implicitly. 

Both symmetries imply that many of the Fourier coefficients vanish in the 
harmonic expansion ( ^ ) of the space-time variances. A complete overview 
IS given m Table |l|. We discuss now the implications of these symmetries 
for a harmonic analysis of HBT radius parameters whose ^-dependences 
are characterised in terms of harmonic coefficients 

^kJ = ^j 4cos(m$)d$, (49) 
Rtj,J = ^ y Rl sin(m$) d<^ . (50) 



2.3.4. Explicit azimuthal dependences at mid-rapidity 

As long as ^-dependent position-momentum correlations in the source are 
weak compared to the global essentially geometrical asymmetry of the fi- 
nite impact parameter collision, the implicit ^-dependence can be neglected 
compared to the explicit one. In this case, all space-time variances ( |4^ ) are 
characterised by their zeroth Fourier components. Inserting the zeroth com- 
ponents listed in Table into the expressions for the HBT radius parameters 
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Table 1. Space-time variances for central and non-central collisions 
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The non- vanishing (x) and vanishing (0) space-time variances for central (6 = 0) and 
non-central (5 ^ 0) collision systems at mid-rapidity (Y = 0) and away from mid- 
rapidity. For non-central collisions, the harmonic coefficients (. . • )5i and {■ .•}^ are 
listed in the three rightmost columns. Away from mid-rapidity, no symmetry ensures 
a different behaviour of odd and even harmonics. 



Hi-- 



-@, one findsBS: 

Rl{K^,<P,Y) 
RUK^,^,Y) 

RUK±,-f,Y) 



m + \ m 

1(052) -(^2)) sin 2$, 

(xz) COS $ , 
— Ixz) sin$ . 



)) cos 2$ , 
)) cos 2$ 



(51) 
,(52) 
(53) 
(54) 
(55) 
(56) 



Two ^-dependent geometrical informations arc contained in these equa- 
tions: 
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(1) Eccentricity of the emission ellipsoid in the transverse planeS 

The HBT radius parameters i?^, i?^ and R^^ are sensitive to the exten- 
sion of the emission region in the transverse plane. They show second 
harmonic oscillations of the same strength ^ ((i^) — (y^)) and thus 
measure the eccentricity of the emission ellipsoid in the transverse 
plane. The second harmonic coefficients of HBT radius parameters thus 
satisfy the rules 

^0,2^ = 0, Rl/ = Q, (57) 

While (|5^) must be always satisfied at mid-rapidity because of the sym- 
metry requirements (see also next section) , deviations from (^^ are an 
unambiguous sign for implicit ^-dependences of space-time variances, 
i.e. for azimuthally dependent position-momentum correlations in the 
source. 

(2) Tilt of the emission ellipsoid in the reaction planeEJ 

The HBT radius parameters R^i and i?^; are sensitive to the off- 
diagonal component (xz). They measure the angle Q by which the 
major longitudinal axis of the emission ellipsoid is tilted away from the 
beam direction, see Fig. ||, 

e= itan"! f-I^^^V (59) 

The tilt is measured from the first harmonic oscillations at mid- 
rapidity. Any deviation from the rule 

Kis=-Ki^- (60) 

is an unambiguous sign for strong azimuthally dependent position- 
momentum correlations in the source. Again, symmetry constraints 
generally require 

Rki=^, (61) 
Ki,i^ = ^- (62) 



2.3.5. Implicit azimuthal dependence at mid-rapidity 

To access the effect of an implicit ^-dependence of space-time variances, 
we expand the radius parameters (^5|)-(|4^) in power of the harmonic 
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reaction plane 



Fig. 5. In general, the main axis of the emission eUipsoid is tilted for finite impact 
parameter collisions by an angle B with respect to the beam axis. This tilt can be 
measured via the first harmonics of the off-diagonal HBT radius parameters R^^ and R^^ 
at mid-rapidity. 



coefRcients The following discussion includes all harmonic contribu- 
tions allowed by the symmetries. It may simplify further if the emission 
function is sufficiently smoothly shaped for higher harmonic coefficients to 
be negligible, 



(63) 



Using the Fourier expansion (^2|) of space-time variances at mid-rapidity, 
one finds: 

1. The zeroth moments of the HBT radius parameters: 



1 



1 



1 



1 



1 7 



(64) 



-P^ {{xi)l + {yi)l) 



Rio = (^'>o ■ 



(65) 
(66) 



These expressions show that contributions to the absolute size of az- 
imuthally averaged HBT radius parameters stem not only from the az- 
imuthal averages (i^ij/)o, but also from the harmonic oscillations of these 
space-time variances. Nevertheless, only those Fourier components of the 
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HBT radius parameters are non-zero, which were found to be non-zero in 
the case without imphcit ^-dependence. This continues to be the case for 

2. The first moments of the HBT radius parameters: 
The only non- vanishing first moments are 

Rii/ = {x~z)o + \ + - pi- mi , (67) 

i?S;,i^ = -(5;5)o + i((iz)^ + (yz)D ■ (68) 

For the case of a static source, these quantities determine the tilt angle ( ^9| ) 
and satisfy i?^; ^ = —R^i ^ ■ According to eqs. ( |67| ) and (|68|), a non-zero 
value of 

Ki.x + Ki.x = mi + mi - /3± m\ (69) 

is an unambiguous sign of dynamically generated ^-dependent correlations 
in the source. 

3. The second moments of the HBT radius parameters: 
The only non- vanishing second moments are 

i?:/ = \ {{fh - i^'h) + \ mi + {i')i) 

- \ {{-X {v% \ mi ' (70) 

Ro/--\mo-mo)+\mi+mi) 

- 13^ {{S:i)l + {yi)l) + i {{iX - {fr:^ + \ ml , (71) 

f mi + imd - i {{-X - (yr^ + \ mi , m 

Ri/-{i')l- (73) 

The equations remain true in the presence of ^-dependent position- 
momentum gradients in the source, while deviations from the ^ = 
~-^s,2^ = ~R'os,2 (|5^ ) are an unambiguous sign for an implicit 
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dependence, namely: 



2 

-(3^ {{xiy, - {yi)l + + {yi)l) , (74) 

I {{^^2 + {yr2)+pi{p)2 - ^ mi - mi) + {iy)i 



- 1(3^ mi + md + \ {{iX mi) > (^s) 

- Ksa = \ m)2 + {f)2) - m)i - mi - (5^)3 - {yi)l) 
~\[{iX~{yX)~{iv)l- (76) 

If the fourth order moments (i^)^ — ijj^Yi {^vYi neghgible and if 
the transverse pair momentum K±^ is sufficiently small, then 

DC 2 DS 2 ^ DC 2 I DS 2 

-^^5,2 - Ros,2 ~ Ro,2 + Ros,2 ■ ['') 

4. The n-th moments of the HBT radius parameters (n > 2): 
In general, the properties of the sin and cos functions imply that the n-th 
order harmonics Rfj Rfj eqs. (|49|), ( |50| ) are built up of m-th order 
harmonics (i^y!^)^ with n — 2 < m < n + 2. This limits the 

number of terms appearing in n-th order expressions. Since it is an open 
question whether the approximation ( |63| ) applies for realistic source^, we 
give these expressions here for completeness. 

For the side-long and out-long HBT radius parameters, only odd har- 
monic terms appear. For n = 3, 5, 7, . . . , we have 

Rol,n^ ^ \ m)n-l - {y^Yn-l) " {^Xn 

+\ m^)n+i + (y^>:+i) , (78) 

Rkn = \ mYn-l - {iX.-l) + \ mYn+l + {^^Yn+l) ■ (79) 

The Fourier decomposition of the other HBT radius parameters contains 
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only even harmonic terms. For n = 4, 6, 8, we have 

(80) 

-l{{-rn,.-{yX,.) + li-yrn.2, (82) 

RiJ = ■ (83) 

Whether these higher harmonic coefficients are numerically important re- 
mains to be established experimentally. 

2.3.6. Implicit azimuthal dependence at forward rapidity 

At non-central rapidity y 7^ 0, additional Fourier components of the space- 
time variances can contribute to HBT radius parameters, as is seen from 
Table |l|. The reason is that at 7^ 0, the symmetries (|46|)"(^8|) do not 
hold. It is remarkable that the zeroth and second moments of the out, 
side and out-side radius parameters, given in eqs. (|6^, (|65[), ([70|), (|7l]) and 
(1?^), do not receive extra contributions at forward rapidity. In particular, 
this implies that deviations from the R% = —R^ = R^s 2^ ^^le of 
(|5^ ) remains an unambiguous test for the presence of angular dependent 
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position-momentum correlations at all rapidities. In addition, the out, side 
and out-side radius parameters acquire first harmonic oscillations away from 
mid-rapidity, 

-\{{-rs + ^i-y)l-{f)l)^ (84) 

+ + , (85) 

+/5i(it)o-i/3x((if)§ + (yt)^) 

+ 1{{-X + Hiy)l~mi) ■ (86) 

In the so-called blast wave modelE^l), there is no correlation between the 
transverse position and the time at which particles are emitted. Hence, the 
space-time variances linear in t vanish. Also, the emission duration (P) 
does not depend on the azimuthal angle. In addition, if the source shows a 
sufficiently smooth azimuthal dependence for the third order terms (i^)^, 
(xy)^, {y'^y^ to be neghgible, then 

RIa ~ Rls ~ "^Rlss ■ (87) 

Moreover, it was observed in a class of model studiesS that should 
be much larger than {xy)\ and {y'^)'i, since asymmetries with respect to the 
beam axis will occur predominantly in the direction of the impact param- 
eter. This translates into the rule 

R'o.i' ■ RIs '■ ^os,i^ 3 : 1 : -1. (88) 

A test of (p7|), ( |88| ) allows to establish whether these additional model- 
dependent assumptions are satisfied. 

The three HBT radius parameters R^ , i?^; and i?^; involve longitudinal 
information and depend on the longitudinal velocity (3i . This leads to addi- 
tional contributions away from mid-rapidity. For completeness, we list here 
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the non-vanishing first and second moments: 

Rl,^ = {~zYi-Wi{~zi)'i^PnP)'i: (89) 

Rh^ = (5'>2 - mm + Pl{pr2 , (90) 

+ 1/3; (2/3^ - 2(H)o - (Si)^ - {m) , (91) 

+ \ + ml) - 1 ((it);; + (y?);) , (92) 

+ iA(2(£t)o-(£i)^-(yt)^) , (93) 

Ki/ = \ imi - mi) + \fii imi - {yi)i) 

+ \ {{x~z)l + ml) - ^ {{S:t)l + {yi)l) . (94) 



2.3.7. Reconstruction of the reaction plane 

The above analysis of HBT radius parameters for non-central collisions re- 
quires the measurement of the angle $ and thus assumes knowledge about 
the event-wise orientation of the reaction plane. This orientation is usu- 
ally measured from the azimuthal dependence of single-particle transverse 
momentum spectrail'0 



d-^p pt dpt dy d<p J 

1 dF' N ^ 

= - ——[l + 2Y,VnC0sn{(t>-^R)]. (95) 

27r Pt dpt dy ^ 

However, the orientation of the true reaction plane can only be measured 
with limited accuracy. Since fluctuations in a finite multiplicity environment 
result in azimuthal anisotropies without geometrical origin, this limited 
accuracy arises largely as a consequence of the basic statistical properties 
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of a mesoscopic system, and cannot be reduced by larger event samples 
or refined measurements. The resulting uncertainty has to be corrected for 
if one aims for a geometrical interpretation of the <i>-dependep£ja^of HBT 
radius parameters. Such corrections are discussed in literatures ^Ir^. 



3. Two-particle correlations from model calculations 

The emission function S{x, K) is not determined uniquely by the correlator 
C{q, K). This is a consequence of the on-shell constraint which implies 
that only a specific time-average over the emission functio n, th e so-called 
relative distance distribution Sk{x), is uniquely measurablal23 

Sk{x) = J dtd{x + f3t,t;K), (96) 
f 4 S(X + ^,K) S(X-^,K) , , 

The direct recon stnictiop of 5j<:(|a;|) from experimental data has been pur- 
sued successfullyESIijiiia. However, due to statistical uncertainties the nu- 
merical inversion of (7(9, fC) — 1 = J (^x cos{q ■ x) S K {x) is complicated. In 
practice, it requires additional model assumptions to achieve convergence. 
Thus most data analyses proceed via model studies. Either they start from a 
model parametrisation of the emission function S{x, K), or they start from 
a dynamical calculation of S{x, K) based on a hydrodynamic or particle- 
cascade based simulation. In this section, we review the main features of 
these approaches. 



3.1. Model parametrisations of the emission function 

The main features of the collision region at freeze-out can be characterised 
by its width in the different spatial and temporal extensions, its collec- 
tive dynamical gradients (usually ascribed to a collective flow field u^j^ix) 
which determines the position-momentum correlations in the source) and 
its random dynamical component (usually ascribed to a local temperature 
T). Model parametrisations of S{x,K) implement these main features in 
a (minimal) analytical ansatz for S{x, K). The model parameters are then 
extracted from a fit to one- and two-particle spectra. 

Example of a model emission function: For illustration, consider a source in 
local thermal equilibrium at temperature T whose extension is given by the 
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transverse width R, space-time rapidity width Ary and proper longitudinal 
emission time tq smeared with a width At. The transverse and longitudinal 
expansion of the collision system results in a longitudinally boost-invariant 
flow profile at freeze-out with a transverse component r?t(r) — Vf^ charac- 
terised by the transverse gradient ryj, 

{x) — (cosh 77 cosh T^i , ^ sinh ?7i , ^ sinh ryt , sinh ?7 coshr/t) . 

rapidity 



(98) 



In longitudinal 
iln[(i + z)/(t- 



Srix,p) 



proper time r = 
z)], this source can be written asol 

2Jr + 1 



V = 



(27r)3 



X exp 



27rAT 

2 



cosh(y 



77' 



- rj) exp 

(t - To 



p ■ u{x) — fir 
T 



2(A7/)2 2(Ar)2 J • ^^^^ 

Here, r labels the particle species which are produced in thermal abun- 
dances with chemical potentials fir- The model emission function (^9|) is 
completely specified by the model parameters T , 77/ , i? , Ayy , Ar , tq , /ir- 
This basic model allows for a satisfactory fit to experimental data from 
the CERN SPsS. 

Overview of models and model extensions: In what follows, we review the 
physics arguments which motivated the study of modifications and exten- 
sions of the parameterisation (p9|): 

(1) Varying transverse density and fiow profiles: 

The functional shape of (^) was varied by replacing the Gaussian trans- 
verse density distribution with a box profilfEj or varying the functional 
dependence of the transverse flow profileE^. This gives further support 
to the general statement that HBT radius parameters are mainly sen- 
sitive to the average r.m.s. width of S{x,K). However, details in the 
functional shape of S{x, K) can leave observable traces in the K± de- 
pendence of the HBT radii: in particular, experimental data from the 
SPS favour a transverse box profile over a Gaussian one. 

(2) Surface dominated versus bulk dominated emission: 

Model ( p9[ ) implements bulk emission, i.e., particles decouple at the 
same average proper freeze-out time tq from all spatial positions in the 
source with a probability proportional to the source density. However, 
if reabsorption of particles by the surrounding matter is signiflcant, 
hadronic freeze-out may proceed via surface evaporation. In analytical 
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parametrisations of S{x, K), such"opaque sources" liave been modelled 
via addition of absorption factorscd. The main outcome of these studies 
is that surface-dominated emission can imply (x^) <C (y^) and this 
makes it possible that < at sufficiently large K±. 

(3) Temperature gradients: 

Models which include a spatially varying local temperature attribute 
position-momentum correlations in the source to a combination of two 
different effects: T{x)and u^{x). This typically introduces two addi- 
tional fit parameterB'0 which characterise the longitudinal and trans- 
verse dependence of T(x). It removes to some extent the constraint 
between temperature and transverse flow which can be exploited when 
fitting (|9^) to a combination of one- and two-particle spectra. Such a 
model was fitted successfully to first data from RHIcEl. 

(4) Model emission functions for finite impact parameter collisions: 



The model (99) has been extended to azimuthaJIy asymmetric tr; 



us 



L51 



verse flow and asymmetric transverse Gaussianii or box geometry 
For the first and second harmonics of HBT radius parameters calculated 
in these models, explicit (f>-dependence dominates over the implicit one 
and deviations from the purely geometrical identities (|5^) and ( ^8| ) 
satisfy the dynamical identities ( |7^ ) and (^), respectively. The latter 
point, however, is mainly a consequence of studying a class of mod- 
els for which t-dependent space-time variances are ^-independent and 
azimuthal deformations are essentially elliptic. 

Generic properties: 

The analytical parametrisations reviewed above were instrumental in 
establishing how the main properties of the two-particle correlator translate 
into specific geometric or dynamical features of the collision region. The 
following generic properties emerge: 

(1) Size and transverse momentum slopes of HBT radii: 

Qualitatively, the main information contained in the absolute size and 
M^-dependence of HBT radius parameters can be illustrated in terms 



of pocket formulas derived in a saddjfi point approximation of (99). For 
the side radius parameter one findsE^I for A?7 = oo and At = 

^?(^^)-T-^^' (100) 
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The size of this radius is proportional to the source size, but it is also 
sensitive to the transverse flow strength 77/ of the source. This illustrates 
that HBT radii characterise only that part of a dynamically expanding 
source which can be viewed through a filter of wavelength K. This 
shrinking effect increases for increasing M±_ proportional to the ratio 
v}/T- The Makhlin-Sinyukov formulaE3 for the longitudinal radius, 

Rf^r^o^^ (101) 



shows compared to (IOC) a stronger Mj^-dependence consistent with the 
stronger longitudinal expansion implemented in (^9|) . Its dependence on 
To is a direct consequence of the assumed longitudinal boost-invariance 
and receives corrections for sources of finite longitudinal extension. 
While quantitative corrections to these analytical expressions can be 
significaniHa, these pocket formulas illustrate qualitatively the interplay 
of geometry and dynamics in determining HBT radius parameters. This 
picture is supported by numerous numerical studies. 

(2) The difference Rl- Rl- 

The main interest in this observablecS lies in its sensitivity to the emis- 
sion duration [see also discussion of eq. (120)] 

Rl-Rl^ al (P) . (102) 

Numerical calculations with a Gaussian density profile typically result 
in a small but positive signal for R^ — R^. For steeper transverse prr 
and particle emission at sufficiently large K± , or for opaque source 
models with surface dominated emission, also negative values can be 
found for R^ — R^. Equation (102) ignores the contribution from the 
(xi) correlation term which vanishes in the model ( |9^ ) but is present 
in hydrodynamic models and Monte Carlo event generators, see below. 

(3) Infiuence of resonance decay contributions: 

Pions from resonance decay s hay e a tendency to be emitted at later 
times and larger distancegH3'E3EZI. For models showing bulk emission, 
their effect on the size of HBT radius parameters is however smallEO. 
This is due to a combination of three effects: i) in models of the type 
( |99| ) , the emission region of the heavier resonances is smaller than that 
of direct pions, ii) the large decay widths of the most abundant reso- 
nances like p's and A's and their non-relativistic velocities imply that 
these decays occur within the emission region of the direct pions, iii) 
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resonances with large lifetime (rj, rj') decay so far outside that their 
decay pions interfere with the directly produced ones on a very small 
relative momentum scale (|q| < 1 MeV) only. This produces a peak of 
the correlation function at |q| < 1 MeV which is narrower than the ex- 
perimental resolution and thus leads to an apparently reduced intercept 
A of the correlator C(q, K) without affecting its shape. Only pions from 
to decays stem from a resonance which is neither sufficiently short-lived 
nor sufficiently long-lived and thus can affect the shape of the corre- 
lator. This spoils a naive core-halo interpretation and contributes to 
non-Gaussian deviations of the two-particle correlator^ 



3.2. Hydrodynamic models 

Hydrodynamic behaviour is an idealised but well-defined limiting case of the 
realistic dynamical evolution of the collision region in heavy ion collisions. It 
emerges as the zero mean free path limit of a particle cascade. In this limit, 
matter in the collision region is treated as an ideal, locally thermalised fluid 
whose dynamics is governed by the relativistic hydrodynamic equations. 
Input for simulations: A hydrodynamic model is fully specified by the equa- 
tion of state and the initial conditions. Typically, the parametrisation of the 
latter models the outcome of an initial pre-equilibrium stage with initial en- 
ergy density estimated from the Glauber, approach to entropy and energy 
production in nucleon-nucleon collisiongl3. 

Freeze-out criterion: The freeze-out criterion, according to which the hy- 
drodynamic simulation is terminated, is another important input in hydro- 
dynamic model studies. Usually, the freeze-out criterion is set by a critical 
energy density or temperature. If the criterion is satisfied in a fluid cell, 
the cell is immediately assumed to freeze-out. Local properties of this cell 
are converted into a thermal ideal gas distribution of hadronic resonances 
with temperature and chemical potential set by the local energy and baryon 
density of the simulation. This leads to a sharp freeze-out along a three- 
dimensional hypersurface and specifles the emission function S{x, K) enter- 
ing the calculation of two-particle correlation functions. There are "hybrid 
models" in which the earlier hot stage of the collision is treated hydrody- 
namically but the hadronic phase is modelled with a Monte Carlo event 
generator codet3. An event generator naturally leads to an emission func- 
tion in a flnite four-volume, see next subsection. 

Successes and problems at RHIC: At RHIC, hydrodynamic simulations 
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Table 2. The main hydrodynamic model calculations with published results 
on HBT correlation functions. Codes follow either the full three-dimensional 
expansion or the two-dimensional expansion in the transverse plane (with as- 
sumed boost-invariance in the remaining longitudinal direction). 



Authors 



ref. 



energies 
studied 



HBT data 
compared to 



Dim. 



HYLANDERkj-y 

Rischke, GyulassyQ 
Zschiesche et al. 

Kolb, HeinzEl 
Hirano, Morita et al. 



NA44EI 
RHIC prediction 

NA49, STAR 
STAR, PHENIX0 



SPS 

SPS, RHIC 
SPS, RHIC 
RHIC 
SPS, RHIC 



(3-l-l)-dim 
(2-|-l)-dim 
(2-|-l)-dim 
(2-|-l)-dim 



NA4' 



90.0, STAR00 (3-Hl)-dim 



compare in general well with the hadronic one-particle transverse momen- 
tum spectra up to « 2 GeV. The major success of this approach is the 
prediction of the size of the measured elliptic flow V2 , as well as the correct 
description of its p_L-dependence for identified pion and proton spectra. This 
indicates that the main contribution to elliptic flow originates in the early 
stages of the collision where the system is very dense and the mean free path 
is close to the hydrodynamic limit zero. However, two-particle correlations 
are determined at freeze-out, where the mean free path (or rather the mean 
scattering time) is grown and a hydrodynamic picture becomes question- 
able. This may be one of the reasons why so far hydrodynamic simulations 
have significant problems in calculating two-particle coirelators which are 
at least in qualitative agreement with experimental dataH, see the following 
discussion. 

Generic properties of hydrodynamic simulations for HBT: An overview of 
hydrodynamic model calculations is given in Table This list is limited 
to studies which include beyond the calculation of onc-particle spectra also 
two-particle correlation functions. 

(1) Freeze-out hyper-surface shows strong outward-temporal correlations: 
The freeze-out criterion implemented in hydrodynamic simulations 
amounts to a sudden switch from a zero mean free path to an infinite 
mean free path approximation. This tends to favour sharp geometrical 
correlations along the freeze-out hyper-surface. In comparison to model 
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sources with emisskm from a finite four- volume, the (it) variance is sig- 
nificantly strongerO. Thus, the difference — does not measure a 
lifetime effect only. This consequence of a sharply localised freeze-out 
hyper-surface may be tamed in hybrid models in which hydrodynamic 
evolution is followed by a hadronic rescattering phaseH. 

(2) Large "lifetime" effect and i?^ > i?^: 

Hydrodynamic simulations lead to sources with very large emission du- 
rations. The size of this lifetime signal depends on the equation of state 
(EOS): a softer EOS results in a more delayed pressure build-up and 
a longer lifetimeElO. Irrespective of model details, the resultingvalues 
for Ro/Rs are generically much larger than the measured resultH. This 
is the main problem of hydrodynamic simulations. 

(3) Resonance decays contribute significantly to HBT radii: 

In contrast to model sources of the type (^^, resonance decay con- 
tributions added to the freeze-out of hydrodynamic simnlations were 
reported to increase the size of HBT radii significantlytJ. This may 
be attributed to the different shapes of the freeze-out hyper-surfaces. 
The homogeneity regions for direct pion and resonance emission are 
the same in hydrodynamic simulations whereas the latter are smaller 
in the model ( |99| ) . Thus even short-lived resonance decay contributions 
tend to increase the hydrodynamic pion source. 



3.3. Monte Carlo event generators 

Event generators are widely used to simulate particle production in ultra- 
relativistic heavy ion collisions. In particular, they allow to study how global 
collective dynamical properties emerge in a mesoscopic system from micro- 
scopic (2-to-2 or 2-to-3 body) interactions. In principle, each event genera- 
tor output defines an emission function from which two-particle correlations 
can be calculated. However, since the event generator output is not a wave- 
function with proper quantum-mechanical symmetrisation, an additional 
prescription is needed of how to relate it to the emission function. There 



is an extensive literature on the conceptual problemli23. In practice, the 
afterburner program of Scott PrattL2l is most frequently used. 

So far, there are only very few calculations of HBT correlation functions 
from event generators, see Table ^. The main conclusion from these calcula- 
tions is that the late hadronic rescattering phase largely determines the size 
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Table 3. Event generator model calculations with published results on HBT 
correlation functions. 



Code 



ref. 



energies 
studied 



HBT data 
compared to 



RQMlfa 
hydrcEi+URQMDi'' 
Humani 

amtpE''' 

MPC0 



SPS 

SPS, RHIC00 



NASsH 

starH 



AGS, SPS, RHIC E859/866, NA44EI STAR0 
RHIC STArH 
RHIC STAR, PHENIxH 



Code not maintained any more, no more recent s1jt«lies available. 
Hadronic rescattering phase dominates HBT radiiEj. 
This code models final state rescattering only. 

This is a multi-phase transport model which includes initial partonic and 
final state hadronic interactions. 

of HBT radius parameters. Also, in contrast to hydrodynamic models, the 
generation of models with relatively small lifetime, satisfying Ro/Rg ~ 1, 
does not appear to be a fundamental problem. While some simulations find 
a ratio Ro/Rs which for large i^x lies betwEen 1.4 and 2.0, inconsistent with 
experimental dataH, other simulationgl3'E3 are consistent with Ro/Rg ^ 1. 



4. HBT measurements 

Two-particle correlations have been measured at all energies from AGS to 
RHIC. Here we give an overview of the experimental situation. 



4.1. Coulomb final state corrections 

Data on two-particle momentum correlations between identical charged pi- 
ous are usually corrected by the experimentalists for the pairwise Coulomb 
repulsion. These corrections are difficult since Coulomb interaction and 
Bose-Einstein interference effects are of similar size and affect the two- 
particle correlator on similar relative momentum scales. Moreover, the used 
correction methods differ between experiments and sometimes even between 
different publications in one experiment. Differences between the used cor- 
rection techniques can change the resulting size of the HBT radius param- 
eter by more than 1 fm and they may affect the iiT^-slope of HBT radii. 
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Thus Coulomb final state corrections are a major source of systematic un- 
certainty in the space-time analysis of correlation measurements. 

The final state Coulomb interaction between two charged particles is de- 
scribed by the relative Coulomb wave-function of the jiaLticle pair, written 
in terms of the confluent hypergeometric function F ljO, 

%T2{r) - r(l + e-^-" e^-^- F {-zr^- 1; z_) , (103) 
z±^\{qr±q-r) = \qr{l± cos 6). (104) 

Here, r = |r|, g = |q|, and 9 denotes the angle between these vectors. The 
Sommerfeld parameter 77 = a/(vrei/c) depends on the particle mass m and 
the electro-magnetic coupling strength e. We write 

f?± = 77: = ±-; , (105) 

where fj, is the reduced mass and the plus (minus) sign is for pairs of unlike- 
sign (like-sign) particles. If particle pairs are emitted from a static source 
at initial relative distance r with a probability >S'stat(''; K), then the corre- 
sponding correlation is given by an average over the squared wave-function 



(IDS) 



Ciq,K) = J d^SsUr^K) l^^^'l^)!' • (106) 
In the case of identical particles, the two-particle symmetrised version of 



(103) should enter equation ( |l0q ). The iiT-dependence of the pair emission 
probability is often neglected when calculating Coulomb corrections. 

The following Coulomb correction methods are based on this starting 
point: 

(1) Point-like Gamow Correction 



For a point-like source S'stat(^) = (5'^-'(r), the correlator (IOC) is given 
by the Gamow factor G{r]) 



$72^(0) 



(107) 



Early studies constructed the corrected like-sign two-particle correla- 
tion by dividing the measured correlator by this Gamow factor 

^(,-,7) (q, K) = Ci-J (g, K)/G(r;_) . (108) 
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(2) Static Finite Size Correction 



The point-like Gamow correction (107) largely overestimates the real 



effect of Coulomb corrections since particles are emitted in reality with 
finite separation r which leads to ajffiaker Coulomb interaction. An 



improvement advocated repeatedlyEHO'Ll is to calculate the correc- 
tion factor for a finite size static source. Typically, a Gaussian ansatz 



(r) cx exp [-r'^/AR^] is chosen in ( |106[) 



J^corr (9) = / d\ S^Ur) l*™/2'(r)P , (109) 
CLrHg, K) = CLas'(9, K)IFl'^l{q) . (110) 



Instead of the analytical emission function in eq. ( 109 ), the particle- 
emitting source can be characterized in terms of a discrete set of phase- 
space points obtained e.g. from a Monte Carlo simulation. In this case, 
the correlation due to particle symmetrization and final state inter- 
actions is usually calculated with a so-called afterburner routine. The 
most widely used afterburner is Scott Pratt's CRABO (CoRrelation 
After-Burner). 



In practice, the value for the source width R in eq. (109) is determined 
iteratively from the extracted HBT radius parameter in the fitting pro- 
cedure. A finite purity of the sample due to misidentified particle leads 
to an overall correlation strength A < 1. This can be taken into account 
by generalisingZa eq. (|11C|) to 



C(rJ(q, K) = il- A) +ACi--\q, K) F^'^liq) . (Ill) 

(3) Correction of like-sign by unlikc-sign correlations 

Rather than to calculate Coulomb corrections for finite size sources, 
one can make use of the fact that unlike-sign correlations receive no 
contribution from Bosc-Einstcin symmetrisation effects but depend on 
Coulomb correlations of the same magnitude (but opposite sign), 

Ci--^ {q, K) = C(+J (g, K) Cl^-J {q, K) . (112) 

Theoretical support for this procedure comes from the fact, that like- 
sign and unlike-sign Coulomb correlations calculated from ( |l06| ) com- 
pensate largely. For point-like sources, e.g., the product of the Gamow 
factors deviates from unity by less than five percent for relative mo- 
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menta q > 8 



137 



A further improvement over (112) is to take this deviation into 
accountO, 

„(--) , _ Citcas {q, K) cLcas (q, K) (^^A) 

^corr.improvcdig' )- G{f]+) G{f]^) ' ^ ' 

This was shown to work with excehent accuracy for a wide range of 
source parameterslZa. 

The effects of finite momentum resolution reduce both G^as^ (g, K) 
and cineas (q, K). As an unwanted consequence, these effects are am- 
plified in the product defining the corrected like-sign correlation func- 



tions (|112|) and (114). An empirical parametrization which takes into 



account finite momentum resolution is discussed below; see eq. (117). 
(4) Experimental parametrisations of Coulomb corrections 

Unlike-sign correlations Cmoas (q, K) were parametrised by the 



functiont£3 

FiQin.) = 1 + (G(r;+) - l) g-^-^/Q" , (115) 
which depends on 



<zi„v = Vq^ - ■ (116) 

The parameter Qq is extracted from the fit. It quantifies a phenomeno- 
logical finite-size correction for large relative momentum. The function 



F(ginv) approaches the Gamow factor (107) for a point-like source, 
Qo — > oo. 

In order to take into account the imperfect purity o f th e sample and 
the finite experimental momentum resolution, CEREsEll parametrised 
the Coulomb correction by 

C^-J {q,K) = {l-\)+X ^ [WK^ (i^coul(ginv) - 1) + 1] . (117) 

In a Monte Carlo simulation of the final momentum resolution, the 
Coulomb correction function i^coui(<Zinv) was obtained by evaluating 
eq. (|l09| ) and reducing F^^f^{q) accordingly. The same Monte Carlo 
simulation determines wkj^ which accounts for the depletion of the pa- 
rameter A due to finite momentum resolution effects. This is chosen 
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such that \wk±^ gives the "true" corrected intercept parameter, which 
then multiphes the correction factor (Fcoui('?inv) ~ !)■ For perfect mo- 
mentum resolution, wk±_ 1 and Fcoui(9inv) — > ^corr(9inv), and the 



prescription (117) agrees with (|lll| 



The above discussion is mainly for static sources and involves Qinv- 
dependent correction factors only. In principle. Coulomb correction effects 
are different for the different relative momentum components, as seen from 



eq. (109). For dynamically expanding sources, a formalism for the calcu- 



lation of Coulomb corrections existso, but it has not been used so far in 
comparison to data. Only one of the correction methods listed above, eq. 



(112), contains some information about expansion effects since it uses the 



measured unlike-sign correlation as correction factor. 



4.2. Experiments at the Alternating Gradient Synchrotron 
(AGS) 

At the AGS of the Brookhaven National Laboratory (BNL) three series 
of collaborations have measured and published results on Bose-Einstein 
correlations in fixed target experiments with beam energies varying between 
2 AGeN and 11.6 AGeW. Due to lack of statistics, measured correlation 
functions were parametrised often by a 1-dimensional parametrisation 

C(gi„v) = l + Ae-'-vflL, (118) 



where the invariant momentum difference ginv is defined in eq. (116). 



E802/E859/E866/E917 These experiments use a rotating spectrom- 
eter (the "Henry Higgins" Spectrometer) which in the E866 upgrade was 
supplemented by a Forward Spectrometer. The acceptance is at or close 
to mid-rapidity. The E802/E859/E866 Collaboration has pubhshed a sys- 
tematic study of the dependence of HBT radius parameters on t rans verse 
mass M±_ — ^ m?- + If ^ , system size and centrality dependencesllH^. The 
last of the series, the E917 experiment, collected HBT data on beam energy 
dependence in Au-|-Au collisions from 6 tg 10.6 AGeN . At the time of this 
writing, these data are still preliminary!! 
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E895 This e xpe riment uses the EOS time projection chamber inher- 



ited from Bevaladli3. BE correlations were measured in 2, 4, 6, and 8 AGeN 



□ □□ 



Table 4. BE correlation data available from AGS experiments. Data are taken at mid-rapidity, are finite size Coulomb corrected 
and are in Bertsch-Pratt parametrisation, unless noted otherwise. The fourth column indicates the analysis frame and the fifth 
column denotes whether the cross-term was included in the fits (yes), not included (no), checked to be and then not included (0) 
or was not applicable since a simpler parametrisation was used (N/A). 



Au-fAu collisions 


impact energy 


coUab.l'et'- 


K± and binning 


frame 


cross-term 


note 


2 AGeV/ 

v^ = 2.4AGeV 


Esgs"* 


0.1-0.3 GeV/c, 3 bins 


CMS 





at midrapidity CMS=LCMS 


4 AGeV 

v^ = 3.1 AGeV 


Esgs""* 


0.1-0.3 GeV/c, 3 bins 


CMS 





at midrapidity CMS=LCMS 


6 AGeV 

v^ = 3.6 AGeV 


Esgs^i"* 

E917"7 


0.1-0.3 GeV/c, 3 bins 


CMS 




N/A 


at midrapidity CMS=LCMS 
preliminary data of -f?inv 


8 AGeV 

V5 = 4.1 AGeV 


Esgs^i"* 

E917"7 


0.1-0.3 GeV/c, 3 bins 


CMS 




N/A 


at midrapidity CMS=LCMS 
preliminary data of Ri^v 


10.6 AGeV 
Vi = 4.7 AGeV 


E917ii^ 






N/A 


preliminary data of -Rinv 


10.8 AGeV 
V5 = 4.7 AGeV 


E877"° 
E877i°8 


0-0.5 GeV/c, 1 bin 
(pj_) = O.lGeV/c 
1 bin 


beam rapidity 


yes 
N/A 


data at beam rapidity 
prelim. -Rinv from K'^ 


11.6 AGeV 
= 4.9 AGeV 


E802107 


3 bins, (rrtx) = 

0.29, 0.4, 0.54 AGeV/c 


YK frame 


yes 




Si-I-X collisions at 14.6 AGeV 


target 


coUab.ref- 


K± and binning 


frame 


cross-term 


note 


Pb 


E814iU« 


1 bin 




N/A 


-Rinv I beam rap., Gamow 




- E802^"''' 


1 bin 


YK frame 


yes 




Au 


E8021'"' 
E8O2/E859I'''' 


3 for TT^ : 

2 for 7r+ 
1 bin 


YK frame 
CMS 


yes 
no 


prelim. K~^K^, Gamow corr. 



3 



S3 

o 
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Au+Au collisions. Eiiblished data exist for the transverse mass dependences 
of correlation radiitij, and the azimuthal dependences in non-central colli- 
sions for the three lower energiesEi^. Results on average phase-space density 
as a functipiLof K±^ and beam energy were shown at the Quark Matter 2001 



conferenceti£ 



All published data are at mid-rapidity. 



E814/E877 Particle correlations at projectile rapidity were measured 
with the Forward Magnetic Spectrometer of the E814/E877 Collaboration. 
The E814 ^a Si-beam experimenillHl while E877 measured Au-t-Au at 



10.8 AGeVliid. E877 determined the pion phase-space density at freeze-out 



for the latter systenJHU. At the Quark Matter 97 conference, they showed 
a direct fit of a Gaussian core-halo source function to K± and ^/-binned 



correlation function and extracted source radii for pions and for kaons 
A summary of AGS experiments is given in Table ^ 



4.3. Experiments at the CERN Super Proton Synchrotron 
(SPS) 

The CERN SPS was used first to accelerate ^^O and 32s nuclei to 200 AGeV. 
Then it was upgraded to accelerate a 158 AGeW ^^^Pb beam. Recently, the 
CERN SPS delivered Pb-beams at lower energies: 40 and 80 AGeV. The 
step-by-step improvement of analysis tools during the CERN SPS heavy ion 
program is clearly seen in the available data. The correlation measurements 
for oxygen and sulphur beams were parametrised first in terms of ginv only. 
The later three-dimensional fits do not include the cross-term ([l8|). Also, 
the Coulomb repulsion was corrected initially by multiplying with a Gamow 
factor which overestimates the repulsion significantly. Improved Coulomb 
corrections, based on averaging the squared Coulomb wave function over a 
finite source size, were only introduced approximately with the arrival of 
the Pb beam. 

NA35 Originally, this experiment used a large streamer chamber in a 
magnetic field to measure tracks of charged particles and their momenta in 
0+Au collisions. The large volume of the detector allowed for the study of 
three rapidity windows between -2.4 and 1.6 in the CMS of the nucleon- 
nucleon collisionEil. Due to small statistics, the correlation function was 
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Table 5. BE correlations data for central oxygen-induced reactions at 200 AGeV/c. All 
correlation functions are constructed from hadron pairs and are corrected for final state 
interactions by the Gamow factor. The parametrisations used by both experiments are 
listed in the last column. They include an additional factor 1/2 in the exponent of the 
correlation function. 



target coUab, 



ref. 



rapidity frame 



A' 



par am. 



C 
Cu 
Ag 
Au 

Au 



WA8dii2 



NASa 



-l<?/lab<l lab 40 - 200 MeV/c 



0.5<j/iab<3.5 
3 bins 



CMS 50 - 600 MeV/c 



ijinv and 
2d: q_i_,qi 



BP 

no cross-term 



parametrised in qi = q\ong and q± - 
C{q) = 1 + A cxp 



1 9 9 1 



qlRl 



(119) 



In contrast to other experiments, this parametrisation used by NA35 has an 
additional factor 1/2 in the exponent. The NA35 detector does not identify 
pions. Thus NA35 always studied hadron-hadron correlations. Wit h th e O- 
induced collisions a measurement of "single-event interferometry"llij was 
attempted. For measurements with the sulphur beam, the detector was up- 
graded with a time projection chamber (TPC) which was crucial in gaining 



good statistics for correlation analysis. In a comprehensive studytHU of S- 
induced reactions with C, S, Cu, Ag, and Au targets, NA35 measured the ra- 
pidity, i^_L, and multiplicity dependence of Bertsch-Pratt correlation radii. 
0+Au results were reanalysed in this work with better statistics. The K± 
dependence of correlation radii was also reported in a letterliiil. Although 
the existence of a sizeable cross-term was first confirmed by NA35liH3, no 
results with the cross-term were published. 



NA49 For the lead beam, this collaboration equipped the NA35 detec- 
tor with four large TPCs which allow for precise tracking of the secondaries 
in the rapidity region 2 < y < 5.5 (values given in the laboratory system 
with ycMS = 2.9 for the 158 AGeV Pb beam). The NA49 detector is able 
to identify particles by a combination of energy loss and time-of-flight mea- 
surements. So far, however, only unidentified hadron-hadron correlations 



□ 



□ 



Table 6. BE correlation data for central sulphur-induced reactions at 200 AGeV. In all cases, final state interactions are 
corrected for by the Gamow factor. None of the parametrisations listed in the sixth column includes the cross term in the BP 
parametrisation. 



target coUab 



ref. 



rapidity 



frame 



par am. 



e 
e 



Al 
Cu 

Ag 
Pb 
Au 



NA35120 
NA44127 

WA80i''2 

NA35120 

NA35120 
NA44127 

NA44126 
NA44124 

NA35120 
WA80"2 



2-5 < yiab < 4.5, 2 bins 

0-5 < 2/lab < 3.5, 3 bins 
yiah ~ 3 

~1 < yiab < 1 

2-5 < j/lab < 4.5, 2 bins 
0.5 < yiab < 4.5, 4 bins 



J/lab ' 
yiab 



0.5 < 



J/lab < 4.5, 4 bins 

-1 < yiab < 1 



CMS 50 - 600 MeV/c BP 

CMS 50 - 600 MeV/c °- BP 

LCMS A'x < 400 MeV/c BP 

lab 40 - 200 MeV/c 2 dim: qj_,qi 

CMS 50 - 600 MeV/c BP 

CMS 50 - 600 MeV/c " BP 

LCMS K± < 400MeV/c 2 dim: g±,qi 

LCMS 2 bins'* BP 

LCMS 2 bins= BP-'' 

CMS 50 - 600 MeV/c "-^ BP 

lab 40 - 200 MeV/c 2 dim: q_i_,qi 



factor 1/2 in the exponent 
factor 1/2 in the exponent 

factor 1/2 in the exponent 
factor 1/2 in the exponent 
factor 1/2 in the exponent 

kaon interferometry 

factor 1/2 in the exponent 
factor 1/2 in the exponent 



" The /fx dependence of correlation radii was measured for 3.5 < yiab < 4.5. 

^ The K± dependence of correlation radii for 2.5 < yiab < 3.5 was also published^^^ . 

The K± dependences of correlation radii were measured in all rapidity bins. 

The two NA44 Kj^ bins correspond to {p±) ~ 150 and 450MeV/c respectively. 

For the BP parametrisation, the two bins correspond to (px) ~ 163 and 246MeV/c respectively. 

Results from other parametrisations are also available'^^''. 
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Table 7. BE correlation data for Pb+Pb collisions at the SPS. Correlation radii were measured as functions of K±^ and the 
cross-term was included in the parametrisations. The fifth column (FSI) specifies how final state interactions are corrected 
for by referring to the corresponding relation in Section 4.1. The column on particle identification (PID) specifies whether 
measurements are for identified pions (YES) or charged hadrons (NO). 



□ 
□ 
□ 



collab.l"ef- 


frame 


rapidity 


par am. 


FSI 


PID 


note 


lvJA44l28,131 


LCMS 


central" 


BP 


(110) 


yes 


also kaon interferometryi3l 


NA45132 


LCMS 


central 


BP 


(117) 


no 


measured for Pb+Au 
@ 40, 80, 158 AGeV/c 


NA49l34,135,137 


FLCMS 


0.1 < ycMS < 2.1 
4 bins 


YKP, BP 


(115)134,135 , 
(112)135, (117)137 


no 


40, 80, 158 ylGeV/ci37 
also kaon 
interferometryi3^ 


WA97143 

□ ° 

WA98"4,Q 


LCMS 
LCMS 


-0.3<3/CMS<0.9 
4 bins 

-0.8<j/CMS<0.2= 


YKP, BP 
BP, YKP 


iterative (110)'' 
iterative (110) 


no 
yes 





" NA44 has a "banana-shaped" acceptance around midrapidity. Small-px data are at slightly forward, high-px data at slightly 
backward rapidity. 

^ In determining the Coulomb correction -Fcorr(<j) from eq. (109) WA97 iteratively used the emission function (99) instead of 
a static Gaussian distribution. 

WA98 has a "banana-shaped" acceptance at slightly backward rapidity. High-px data are further away from midrapidity. 
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Data on the centrality dependence of correlation radii. 



coUab 



ref. 



system and energy 



rapidity 



E802E£3 



NA44 
NA4,i3' 



NA4C 



136 



WA9' 



STARl 



Si+Al @ 14.6 AGeV 
Si+Au @ 14.6 AGeV 
Au+Au @ 11.6 AGeV 



NA35E3 S+Ag @ 200 AGeV 



S+S, S+Ag, S+Pb 
@ 200 AGeV 

Pb+Au @ 40, 80, 158 AGeV 
Pb+Pb @ 40, 80, 158 AGeV 
Pb+Pb @ 158 AGeV 



Au+Au O = 130 AGeV ii 
= 200 AGeVll£^ 



PHENIxtia Au+Au @ = 200 AGcV 



mid-rapidity 

-2.5 < J/CMS < 1-5, 4 bins 
mid-rapidity 

mid-rapidity 



J/CMS = 0| 



H, 1. 



-0.3 < J/CMS < 0.9 

mid-rapidity 
mid-rapidity 



are available. For Pb+Pb collisions at 158 AQeY . the rapidity and 
dependence of YKP parameters was publishedHil. Also, a rather compre- 
hensive compilation of preliminary data at all rapidities and in both YKP 
and Bertsch-Pratt is availablelHa. Preliminary resultal^S't^ exist for h^h^ 
correlations in 40, 80, and 158 AGeV Pb+Pb collisions with different ccn- 



tralities. Results for kaon-kaon correlaii 



collisions were published very recentlyt^a. These latter studies are the only 



pns in central 158 AGeV Pb+Pb 



ones based on "global tracking" for which tracks from all TPCs are matched 
before constructing the correlation function. 



NA44 This experiment is based on a focusing spectrometer. Its narrow 
acceptance is around mid-rapidity and depends in detail on the detector set- 
ting which may be varied. The performance is optimised for particles with 
small momentum difference and good particle identification is achieved, 
the S-beam, one-dimensional correlation functions for S+Pb coUisiond 
and a kaon interference studyt^ were published first. Bertsch-Pratt radii 
without the cross-term were published in 1995I1H3 and a dedicated paper was 
written on their M± dependencellHl. NA44 studied collisions with other 
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collab.'^®^- 


system and energy 


rapidity 


NA44 




S+Pb @ 200 AGeV 


mid-rapidity 


WA9^ 


Pb+Pb @ 158 AGeV 


-0.8 < ycMS < 0.2 








(ycMs> = -0.2 


STAR 




Au+Au @ = ISOAGeV 


mid-rapidity 



Table 10. Data on the azimuthally dependence of HBT radius param- 
eters. All data are taken at central rapidity. 



coUab 



ref. 



i-dependent radii system and energy 



E895 



all six radii {(f>) Au-|-Au @ 2, 4, 6 AGeV 



STARO Ro, Rs, Rls, Rf Au+Au @ ^ = 130, 200 AGeV 



targetslilil: S-fS, S-f Ag, S+Pb. More recently an mwstigation of three- 
particle correlations in S-|-Pb systems was publishedtH; . 



There are_two papers with da ta fr om the 158 AGeV Pb-|-Pb collisions, 
one for pionlH3, and one for kaonli^ interferometry. These study the AI± 
dependence of the BP correlation radii within the limited acceptance of the 
detector. 



NA45-CERES In 1998, the CERES collaboration upgraded their de- 
tector with a time projection chamber with radial drift field. This al- 
lows for interferometric studies. CERES measured correlations of non- 
identified hadfons in Pb+Au collisions at 40, 80 and 158 AGeV with various 



centralities!£2 



WA80 The correlation analysis of WA80 is based on the so-called plas- 
tic ball detector which has coverage in the target rapidity region and at 
low p± < 220MeV/c. For the oxygen beam they used C, Cu, Ag, and 
Au targets. The sulphur beam was collided with Al and Au targets, the 
proton beam at 450 AGeV/c was collided with C and Au targets. An ear- 
lier analysis of target dependence of the observed HBT radii for 0-induced 
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Fig. 6. Summary of the -dependence of correlation radii measure 
SPS in Pb+Pb cpllifiions at-^idrapidity. iBata are taken from NA4' 
preliminary NA49t£3, WAQTtla and WAgif • 



,t the CRB^f 
CERESE3, 



reactionallJ was superseded by a study including also S and p as projectiles 



in which the detector performance correction was better understoodlllj. The 



latter paper also shows a simple model fit to the correlation function with 
Coulomb correction (IfO) instead of Gamow factor multiplication (f08). 
They used an addilimial factor of 1/2 in the correlation function like NA35 
did. Another worktlll analyses these data in the context of intcrmittcncy. 



WA98 For the lead beam runs, this collaboration made use of the 
plastic ball calorimeter. To measure charged particles the detector in- 
cludes a two arm tracking spectrometer with "banana-shaped" tt" ac- 
ceptance around mid-rapidity. WA98 published a study of the K_\^ de- 
pendence of cor relation radii for identified tt" in Ber tsch-Pratt and YKP 



parametrisation 



144, 146 



Three pion interferometryli 



145,146 



and pion phase- 



space densitytfj were also studied. 
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Fig. 

at RHIC in Au+Au collisions. Data are taken from STAR|i£^t£3 and PHEND 
Data for y/s = 200j4GeV are taken from transparencies shown at the Quark Matter 2002 
conference. 



7. Summary of the K i -dependence of correlation ra dii . at m id-rapidity m«as,u¥ed 



WA97 To measure BE correlations, WA97 used a silicon telescope 
which provides precise tracking of produced particles within the mag- 
netic field. They do not identify particle species but sample the correla- 
tion function with h^h~ pairs. The data are mainly presented in the YKP 
parametrisation but consistency checks with the Bertsch-Pratt form were 
performedtlS. Transverse momentum and rapidity dependences of the cor- 
relation radii were investigated in the acceptance window —0.3 < y < 0.9 
(in CMS) and 0.2 GeV/c < pj. < 1.3 GeV/c. 



4.4. Experiments at the Relativistic Heavy Ion Collider 
(RHIC) 

Three of the four collaborations at RHIC published results on BE inter- 
ferometry: STAR, PHENIX, and PHOBOS. Data were taken from Au+Au 
collisions at CMS energies of 130 and 200 AGeV. 
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STAR uses a time projectioiLchaniber in a solenoidal magnetic field. 
The analysis for ^/s — 130 ^GeVt^ covers in three bins the if^-dependence 
from 0.125 GeV/c to 0.45 GeV/c within a rapidity region |?/cms| < 0.5. 
Preliminary data, at ^/s = 200 AGeV range in four bins up to K± k, 
0.52 GeV/ai£ilIl|e centrality dependence of the HBT radii at both energies 
was studiedt^ll^. The azimuthal dependence of the ifj^-integrated HBT 
radii ff^ . i?^ and R^^ was also analysed at both energies in minimum bias 
evental^j. So far, STAR determines the orientation of the reaction plane 
but it does not determine the direction of the impact parameter. Thus, the 
azimuthal dependence of i?^; and i?^; cannot be measured. Finally, a first 
study of three pion correlations indicates that the source is fully chaotic 
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PHENIX extends the results of STAR in a pseudo-rapidity window 
\'q\ < 0.35 to higher K^- 0.2 < Kj_ < 1.2 GeV/c at the lower and 
0.2 < K± < 2 GeV/c at the higher CMS energy. The particle momentum 
is measured by a drift chamber and a pad chamber. At = 130v4GeV, 
correlation radii for b oth positive and negative identified pions are deter- 



mined in three K± binslilH. For collisions at y/s ~ 20 Q AG eV. a much better 



statistics allowed to split the pairs into nine K± binsElj. The centrality de- 
pendence of the correlation radii was also studied. Correlation radii from 
kaon-kaon correlations did not showllla a simple M±_ scaling with the tttt 
radii in contrast to expectations from certain hydrodynamically motivated 
parametrisations of the freeze-out state of the fireballEj. 

PHOBOS presented so far two sets of BP correlation radii for the 15% 
most central -^/s — 200 AGeV Au-|-Au collisions. Data are for one bin 
from 0.15 to 0.35 GeV/c and for 0.2 < V < 1.5. One set was measured with 
7r+ pairs, the other one with vr^ pairglla. 



4.5. Discussion of the data 

4.5.1. Size and transverse momentum dependence of HBT radii 

The out-, side-, and longitudinal HBT radius parameters vary typically 
around 5 — 6 fm at small transverse pair momentum K± and decrease with 
increasing K^^. Their absolute size shows no significant dependence on beam 
energy (see Figs. ^ and For data from the CERN SPS, the Xj^-slope 
of the longitudinal radius parameter is steeper (see Fig. ^ while all radius 



HBT from AGS to RHIC 



47 




Fig. 8. The awd-iKx dependence of Ro and-Ss . Data without error bars ar 
marised from ESgslliJ (Au+Au), NA45-CERESE£3 (Pb+Au), and STAR 7r+7r+E£^ 
( Au+ Au) . STAR results for ^ = 200AGeV are taken from transparencies shown at the 
Quark Matter 2002 conference. 
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Fig. ft^ — I The ^/s and ii" jr-d«pendence of i?;. Da^tt- 



E895E. 



, NA45-CER,Egi£3, and STAR ■k+tt+ 



out error bars are summarised from 
STAR results for ^ = 200^GeV 



are taken from transparencies shown at the Quark Matter 2002 conference. 



parameters measured at RHIC show approximately the same slope (see 
Fig. 0). 

As explained in section 2.2.2, the HBT radius parameters of an ex- 
panding source correspond to the width of the X^-dependent region of 
homogeneity. This is smaller than the width of the entire collision region. 
The -ftr_L-slope of the HBT radii is a measure of the collective dynamical 
expansion. This picture can be illustrated by the pocket formulas (100), 
(101) and is supported by many model comparisons. If SPSdata are fitted 
by a model with Gaussian transverse density distributionBH, this leads 
to a radius of the entire collision region i? « 7 fm. To put this number 
into perspective, we relate the two-dimensional rms width of the collision 
region, ^j^J^Jlf'^'^ = \/ (i^ -I- y^) — \plR ~ 10 fm, to the two-dimensional rms 
widths of a cold lead nucleus. The hard sphere radius R^^ = 1.2 A^^^ fm is 
for lead = 7.1 fm, and the corresponding two-dimensional transverse 
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rms width is r^°^f^^ = y^{x^ +y^)ph = ^iri ^''^ '^^^^ 

favour a model with a transverse box density profiletJ over a Gaussian den- 
sity profile. The box radius is f — 12 fm. Irrespective of the transverse 
profile, one concludes that during the collision the system has expanded by 
a factor w 2 from the transverse size of the overlapping cold lead nuclei to 
the transverse extension at freeze-out. 

All available data are subject to significant systematic uncertainties. 
Additional uncertainties arise when comparing data from fixed target and 
collider experiments. In view of the rather mild changes of HBT radii be- 
tween SPS and RHIC, this makes it difficult to assess to what extent the 
dynamical interpretation given above changes from SPS to RHIC. A first 
analysis of RHIC dataE3 argues in favour of a more extended source with 
larger transverse flow, thus supporting the picture of a more vigorous trans- 
verse expansion at higher centre of mass energies. 



4.5.2. Ro/Rs 

The main interest in the quotient or difference of the two transverse HBT 
radius parameters ( [l5| ) and (|6|) lies in a model-deapmident argument, that 
the emission duration (P) can be extracted from Be3 

RliK) ~ RliK) ^ (3l{P) . (120) 

This statement is based on two model-dependent assumptions. First, the 



term —2Pj_{xt) should be negligible compared to (120). This assumption 



however, can be violated in models with strong expansion. Second, the 



difference (i^) — (y^) should be negligible compared to (120). This latter 
assumption can be violated at sufficiently high K± , in particular in models 
for which particle emission peaks close to the surface due to dynamical or 
opacity effects. 

Many model calculations predict Ro/Rg 3> 1. In particular, Rischke and 
GyulassyO emphasised that this would be an unambiguous signal of an 
equation of state which is sufficiently soft in the phase transition region to 
result in a significantly delayed build-up of transverse expansion. This would 
result in a large lifetime effect, Ro/Rg ^ 1-5. In contrast, data indicate 



values Ro/Rg < 1.1 even at RHIC, see Fig. As mentioned in section 3.1 



one can think of physics effects which result in Ro/Rg < 1. This issue is 
59] 



presently under studyE 
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0, 



Fig. 10. Ther-y^ and K^^ depend nnem of i Ro/ Rs ■ Data are summarised from E895I 
NA45-CEREslii3, and STAR vr+Tr+ESEB. STAR results from ^ = 200ylGeV are taken 
from transparencies shown at the Quark Matter 2002 conference; their error bars may 
be underestimated. 



4.5.3. Average phase-space density 

Idea and formalism: The spatial average of the phase-space density of pions 
at times later than freeze-out tf 



Jd^xp{x,p,t> tf) 
J d3xf{x,p,t> tf) ' 



(121) 



can be exp ressed in terms of the observable one- and two-particle 



spectr 



{mK^,Y) = 



V,siK^,Y) = 



dN 



1 1 

TtEdY M^dMi_ d$ V,s{Kj^,Y) 
j d^q{C{q,K)-l) . 



(122) 
(123) 
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E895 Negative Pion Phase Space Density 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Pj (GeV/c) 



Fig. 11. Preliminary data 
on the average phase-space 
density as measured by the 
E895 Collaboration for dif- 
ferent AGS projectile en- 
crgftesi as a function of 
Pj MH . The lines show 
Bose-Einstein distributions 
of given temperature. This 
does not account for the ef- 
fects of expansion. 



The phase-space density is determined by the average number of pions with 
given momentum (the non-invariant spectrum) divided by the volume in 
which they are contained. The factor A~^/^ corrects for the "purity" of the 
sample: it ensures that only directly produced pions and not those coming 
from resonance decays are taken into account (see discussion at the end of 



section 2.1). As long as lifetime effects are small (which is consistent with 
all data measured so far, see Sec. 4.5.2 ), the volume is given by (123). For 
the Cartesian BP parametrisation, it takes the form 



(124) 



Thermal model as a reference: For reference, it is customary to compare 
the measured (/) to that of a thermal distribution of a given temperature 
and chemical potential. Since (/) is a spatial average, it averages for an ex- 
panding source over homogeneity regions which move relatively to the point 
of maximum emissivity. This leads to deviations from a naive Bose-Einstein 
phase-space density which for the case of a boost-invariant longitudinally 
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Fig. Ip^^ The average pion phase-space density measured by NA44@, NA49Ei, and 



O 



WA9I 



in 158 AGeV Pb+Pb collisions at mid-rapidity. 



expanding pion source with transverse expansion can be written asE^I 

{f)iP±) = {T.Mp±)^ j (E("-1)^"(^^^)) ' (125) 

^0 

xIo(np^sinhC(r)/T) Ki(nmj^ coshC(r)/T) . (126) 

Here, the transverse rapidity profile is characterised by C^(r) and the trans- 
verse density profile by ^{r)/T. For illustration, the transverse geometry 
can be modelled e.g. by a box profile 



[box] 

or a Gaussian profile 



/is for r < i?box 
—oo otherwise' 



(127) 



[Gauss] /i(r) = — T 



^^Gauss 



(128) 
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Fig. 13. The average pion phase-space density as measured by STAR CoUaboratiorJi^ 
for Au+Au coUisions at = 130AGeV/c. Different symbols correspond to centralities 
(from top to bottom) 0-5%, 5-10%, 10-20%^0-30%, 30-40%, 40-50%, and 50-80% of the 
total cross section. Curves are fits to eq. (125). The same thermpii model also reproduces 
the single-particle spectrum. Triangles represent the NA49 dataO for 158 AGeV Pb-|-Pb 
collisions at mid-rapidity. 



The chemical potential /i(r) introduced here are position-dependent and re- 
sult in a non-uniform density profile; /is or fiQ are the values in the centre 
of the fireball. A certain spatiaL average of this chemical potential deter- 
mines the particle multiplicitiea23. For comparison of the present formalism 
to AGS data, the assumption of boost-invariance entering (125) has to be 
modified by a rapidity cut-off. 

Data: The average phase-space density was studied first by the E877 experi- 



ment for Au-|-Au collisions at the AGSli23'Ei£l in the projectile fragmentation 



region where it was found to decrease with increasing rapidity. The E895 
collaboration measured the average phase-space density at mid-rapidity at 



2, 4, 6, and 8 AGeV and observed its increase with the collision energylH^ 
(Fig. 0). 

At SPS energies there is an extensive compilation of phase-space den- 
sities for various collision systems based on data of the NA35/NA49 
collaboration^, the NA44 measurements of S-l-S, S-|-Pb, and Pb-|-Pb 
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Fig. 14. Dependfftief; of the 'jfr«eze-out volmtp" V/ = (27r)^/^_Rg_R( on Data are 
measured by ESQaliJ, CERESEB and STARE53 for low-px pions at mid-rapidity. 



collisiondlHd, and the Pb+Pb data measured by WA98. For Pb+Pb colh- 
sions at 158 AGeV, the results of these experiments are mutually consistent, 
see Fig. 

At RHIC energy, preliminary data indicate that the average pion phase- 
space density is significantly higher than at SPS. Moreover, the STAR 
collaboratimi reported a strong dependence of the phase-space density on 
centralityE-^ (Fig. ^3|). The origin of this increase, as well as the apparently 
different p^-dependence of (/) at RHIC and SPS, is currently under study. 



4.5.4. Energy and multiplicity dependence 

Despite the weak dependence of HBT radius parameters on centre of mass 
energy, the volume Vj from which particles decouple shows an interesting 
non-monotonous behaviour, see Fig. [l^. Here the volume is estimated as 
Vf DC R^Ri. The radius Ro is not used in calculating Vf because it contains 
contributions from the temporal extent of the source, see eq. (|l6|). The 
estimate Vf cx R^Ri does not take into account the effects of the expansion 
on the HBT radii. Figure |lj shows that the volume decreases gradually 
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within the AGS energy range, reaches a minimum around the highest AGS 
energies and then increases monotonously by almost a factor 2 up to the 
highest RHIC energy. 

At approximately fixed centre of mass energy at the SPS, the same 
freeze-out volume was found previously to grow linearly with the charged 
particle multiplicity per unit rapidity, see Fig. Consistent with this find- 
ing is the centrality dependence of the freeze-out volurne at SPS energies 
which grows linearly with the number of participantsll23. These two ob- 
servations support the conjecture that freeze-out occurs at a fixed particle 
density. This is, however, contradicted by the non-monotonous energy de- 
pendence of Vf. A linear relation between freeze-out volume and particle 
multiplicity does not hold. 

Particle density, and thus particle multiplicity is certainly important in 
characterising the freeze-out condition, since it affects the hadronic escape 
probability from the medium. However, chemical composition, collective ex- 
pansion and the momentum of thpescaping particle are other factors which 
determine this escape probability^. To illustrate this, one can consider e.g. 
the mean free path of a pion at freeze-outi 



(Jl33 



'^mfn = 2rip+pCr7rAr + 3n^-cr^^ = 2 ct^at + 3— ^CTtttt ■ (129) 

Vf Vf 

While the freeze-out volume Vf and the numbers Ni of particles of species 
i contained in Vf both depend significantly on the centre of m ass energy, 
the mean free path ( 129 ) is approximately y^-independent 



L33 



To understand in more detail how the interplay of different properties of 
the collision region determines the freeze-out volume, a realistic freeze-out 
criterion is required. A good starting point is the particle escape probabilitv 
from the hot and dense but rapidly expanding collision regiono'E§ 

n..P,r)^e.p(-f + , (130, 

Here, v is the velocity of the escaping particle and Tl{x,p) denotes the 
scattering rate which is defined as the inverse of the mean time between 
collisions for a particle at position x with momentum p. Freeze-out at differ- 
ent centre of mass energies is then assumed to occur when the probability 
V reaches a characteristic value. 

Recently, the scattering rate TZ{x,p) was calculated for a full hadron 
resonance gas with chemical composition corresponding to SPS and RHIC 
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Fig. 15. The "freeze-out volume" (Vj oc \/?:B?^Ri) as a functionr-e^ the number of 
produced negative particles per unit of rapidity. Data are by NA35E2a for different S- 
induced reactions at 200 AGeV beam energy. 

energiesiil. The main observation is that the scattering rate shows a sig- 
nificant momentum dependence suggesting that particles of different mo- 
menta are emitted at different times. This effect is neglected in hydiody- 
namic simulations which are based on the Cooper-Frye prescriptionca for 
freeze-out along a sharp three-dimensional hypersurface. Deviations from 
this Cooper-Frye prescription, i.e. freeze-out along finite four volurn£s may 
affect the transverse momentum slope of HBT radius parametersEJ, and 
the momentum dependence of the freeze-out volume. The role of a change 
in the chemical composition from SPS to RHIC was found to be relatively 
small in spite of the large increase of pion phase-space density (Fig. p^ ). 
This is a consequence of the small pion contribution to the total scatter- 
ing rate, resulting from the comparatively small pion-pion cross-section. In 
contrast, collective transverse expansion gradients affect the freeze-out vol- 
ume significantly. The reason is that an increase in the scattering rate at 
freeze-out can be compensated by stronger transverse flow gradients which 
lead to a faster density decrease in the collision region thus keeping the 



opacity integral in the exponent of (13C) constant. The possible effect of 
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Fig. 16. Azimwthal dependence of HBT radius parameters as published by the E895 
Collaboration for Au+Au collisions at 2 AGeV. Curves correspond to a static source 
((P) = 0) according to the equations (|5l[)-(|56[). 



the flow gradients on the -yi-dependence of the freeze-out volume (Fig. 
remains to be studied. 



4.5.5. Azimuthal dependence of HBT radius parameters 

Two years ago, the first measurements of the ^-dependence of HBT radius 
parameters were pubhshed by the E895 CoUaborationliiJ for beam ener- 
gies of 2, 4 and 6 AGeV in semi-peripheral Au-|-Au collisions at the AGS. 
Results at all three energies show a sizeable first order harmonics in the 
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^-dependence of i?^; and i?^;, and comparatively small second order har- 
monics in the fully transverse HBT radii i?^, i?f and i?^^. The longitudinal 
radius parameter and the intercept A are consistent with a ^-independent 
ansatz. According to (|59|), the first order harmonics allow to reconstruct the 
angle 8 by which the emission ellipsoid is tilted out of the beam axis, see 
Fig. 1^. At lower AGS energies, this angle is with O « 30° surprisindv large. 
This value is consistent with RQMD transport model simulationsEZl. Inter- 
estingly, the spatial tilt is found to point in the direction opposite to the 
directed flow in momentum space. This indicates that at lower AGS ener- 
gies pion reflection from the bulk of the matter rather than pion absorption 
by this matter is at the root of the observed direct flow signal. 

With increasing centre of mass energy, a longitudinally approximately 
boost-invariant region develops around mid-rapidity. As a consequence, the 
tilt angle O of the emission ellipsoid is expected to decrease with increasing 
y/s. However, there is so far no measurement of and i?^; at higher 
energies, which would be required to establish this effect experimentally. 

The ^-dependence of the fully transverse radius parameters i?^, R^^ 
and i?g is easier to measure than that of R^i and i?^; : while the former re- 
quire the event-wise reconstruction of the orientation of the reaction plane, 
the latter require in addition the direction in which the impact parame- 
ter points. Due to this complication, at RHIC first preliminary data are 
available for the fully transverse radius parameters onlyll£3. These data are 
expected to contain information about whether the spatial orientation of 
the source at freeze-out is in-plane or out-of-plane. However, statistical and 
systematic uncertainties in these preliminary data are still too large to draw 
conclusions. A significant improvement in statistics is expected within the 
next run. 
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